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Long-Term  Goals 

Ocean  waves  are  an  important  aspect  of  upper  ocean  dynamics,  in  particular  on  shallow  continental 
shelves  and  in  coastal  areas.  The  long-term  objective  of  this  work  is  to  advance  modeling  capability  in 
such  coastal  areas  by  improving  model  representations  of  effects  associated  with  nonlinearity,  medium 
variations,  and  dissipation. 

Objectives 

The  specific  objectives  of  the  present  work  are  1)  to  develop  and  implement  an  efficient  and  scalable 
approximation  for  the  nonlinear  quadmplet  source  term,  2)  to  develop  a  generalized  nonlinear  source  term 
that  is  accurate  in  water  of  arbitrary  depth,  3)  to  develop  improved  nonlinear  closure  approximation  for 
triad  nonlinear  interactions  in  shallow  water,  4)  to  develop  and  implement  a  new  stochastic  model  for 
coherent  wave  interference  induced  by  varying  depth  and  currents,  and  5)  improve  representations  of 
dissipation  by  wave  breaking  and  wave-bottom  interactions  in  shoaling  waves. 
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Approach 

Modern,  third-generation  (3G)  wave  models  are  based  on  an  action  balance  (or  radiative  transfer) 
equation,  which  describes  the  transport  of  wave  energy  (or  action)  through  a  slowly  varying  medium  and 
time.  In  Lagrangian  form  (for  convenience)  this  balance  equation  can  be  written  as 

- =  Sin  (k)  +  Sds  (k)  +  Ssc  (k)  +  Snl  (k)  (!) 

at 

where  N(k)  is  the  wave  action  at  wavenumber  vector  k  and  t  is  time.  The  forcing  terms  on  the  right- 

hand  side  are  referred  to  as  source  terms  and  account  for  the  input  of  energy  by  the  wind  (S;n),  spectral 
redistribution  of  energy  through  scattering  by  seafloor  topography  (Ssc)  or  through  nonlinear  wave-wave 
interactions  (Sn[) ,  and  dissipation  of  wave  energy  (Sds)  through  e.g.  breaking  or  bottom  friction. 

In  this  study  we  will  develop  and  improve  the  source  terms  for  nonlinear  interactions  Sni  and  energy 
dissipation  Sds,  to  account  for  effects  of  finite  depth  and  shallow  water,  and  to  ensure  a  consistent  and 
smooth  model  representation  of  wave  evolution  from  deep  to  shallow  water. 

Nonlinearity 

We  will  develop  theory  and  modeling  capability  to  transport  (nonlinear)  cross-correlators  through  regions 
of  variable  depth.  To  allow  modeling  of  wave  propagation  from  deep  to  shallow  water,  we  will  modify 
the  nonlinear  source  term  to  account  for  changes  in  relative  water  depth  [Janssen  et  al.  2006],  and  develop 
an  improved  closure  approximation  for  nearshore  wave  propagation  [Janssen,  2006],  To  allow  for  more 
efficient  evaluation  of  nonlinear  four-wave  interactions  in  operational  wave  models,  we  develop  a 
Lumped  Quadruplet  Approximation  (LQA)  by  rewriting  the  exact  WRT  method  to  reduce  the  number  of 
effective  configurations  by  lumping  contributions  along  the  loci-integration. 

Dissipation 

We  test  wave  dissipation  in  muddy  environments  by  comparing  model  simulations  with  an  operational 
wave  model  to  recent  field  observations  of  waves  propagating  across  a  muddy  seafloor  for  a  range  of 
wind  and  wave  conditions. 

Field  data 

We  will  analyze,  prepare,  and  disseminate  selected  field  experimental  data  sets,  collected  by  the  Pi’s  and 
colleagues,  to  the  project  teams  for  the  purpose  of  validation  and  calibration  of  new  model  developments. 
The  data  will  be  made  available  through  a  dedicated  website. 
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Work  Completed 


Development  of  a  Lumped  Quadruplet  Approximation  (LQA) 

A  first  step  in  the  inter-comparison  of  exact  methods  has  been  carried  out  using  the  WRT  method  as 
implemented  by  Van  Vledder  (2006)  and  the  RIAM  method  of  Masuda  (1980).  The  results  of  the  inter¬ 
comparison  were  presented  at  the  13th  International  Workshop  on  Wave  Forecasting  and  Hindcasting, 
Banff,  Canada,  (Van  Vledder  and  Hashimoto,  2013).  Results  of  Lavrenov’s  Gaussian  Quadrature  Method 
(see  Lavrenov,  2001),  as  implemented  by  Gagnaire-Renou  et  al.  (2010)  are  now  being  processed  and 
added  to  the  inter-comparison  study.  Moving  forward,  the  recently  introduced  exact  method  of  Onorato 
and  Janssen  (2013)  and  analytic  solutions  (Dungey  and  Hui,  1979)  for  narrow-peaked  spectra,  will  be 
added  to  the  comparison. 

The  LQA  and  the  Generalized  Multiple  DIA  have  been  implemented  in  subroutine  form  and  included  in 
the  research  wave  model  WAVETIME  of  Van  Vledder.  Since  each  approximation  of  the  nonlinear  source 
term  will  affect  other  source  terms  as  well,  the  focus  in  the  assessment  of  the  new  approximation  is  on  the 
reproduction  of  integral  model  results,  rather  than  the  ability  to  represent  the  details  of  the  nonlinear 
transfer  rate  for  an  arbitrary  spectrum.  The  LQA  and  the  GMD  are  now  being  implemented  in  the  SWAN 
model  to  evaluate  its  performance  in  growth  curves  analysis  and  field  cases. 

Analysis  and  dissemination  field  observations. 

We  have  developed  an  interactive  website  for  the  dissemination  of  field  data  sets  to  the  NOPP  teams  (and 
the  community  at  large).  The  website  (www.waveserver.org)  includes  an  interactive  Google  Earth 
implementation  that  allows  the  user  to  visually  browse  the  various  datasets  (see  Figure  1).  Thus  far  we 
have  made  the  following  data  sets  available  through  this  website: 

•  1999  ONR  Shoaling  Waves  Experiment  (SHOWEX),  wave  transformation  across  a  wide 
continental  shelf; 

•  2003  ONR  Nearshore  Canyon  Experiment  (NCEX),  refraction  over  nearshore  bathymetry; 

•  2004  SAX04,  Florida  Gulf  Coast  (ONR  Ripples  DRI),  wave  attenuation  due  to  bottom  friction; 

•  2007  Martha’s  Vineyard  Experiment  (ONR  Ripples  DRI),  nearshore  wave  propagation  over 
inhomogeneous  bottom  decomposition; 

•  2008  ONR  Waves-over-Mud  MURI  (Louisiana),  wave  propagation  over  muddy  seafloor; 

•  2009-2010  High-resolution  Air-Sea  experiment,  drifter  observations  of  coastal  wave  evolution; 

Of  particular  interest  to  the  NOPP  team  is  the  recently  added  SHOWEX  data,  which  captures  wave 
evolution  across  the  wide  shelf  of  the  Mid- Atlantic  Bight,  in  the  region  between  Cape  Hatteras  and  the 
entrance  to  the  Chesapeake  Bay,  and  includes  detailed  measurements  of  hurricane  swells  (Ardhuin  et  al., 
2003a, b)  and  fetch-limited  wave  growth  (Ardhuin  et  al.,  2007). 
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Figure  1  Left  panel:  Snapshot  of  www.waveserver.org.  This  website  is  developed  as  part  of  this  NOPP  project  to  disseminate  field 
observations  collected  by  the  Pis  and  colleagues,  with  the  NOPP  teams  and  the  community  at  large.  The  home  page  includes  an 
interactive  Google  Earth  window  that  allows  users  to  browse  through  the  various  datasets.  A  link  to  the  dataset  appears  when  the 
user  clicks  on  the  buoy  icon  in  the  panel.  The  right  panel  shows  a  zoom-in  of  the  Google  Earth  window  with  labels  for  each  of 
the  experiment  sites  for  which  data  is  available. 


Testing  of  benchmark  deterministic  models 

To  develop  a  benchmark  model  for  further  validation  of  the  shallow-water  closure  approximations 
developed  in  this  project,  we  have  implemented  and  tested  a  non-hydrostatic  model  (SWASH)  to  capture 
nonlinear  wave  propagation  in  a  dissipative  surf  zone  (see  Smit  et  al.  2014a).  Non-hydrostatic  models 
provide  an  efficient  means  to  model  nonlinear  wave  propagation  and  the  results  are  very  promising  (see 
Figure  2).  However,  for  the  sake  of  efficiency,  non-hydrostatic  models  assume  a  single-valued  free 
surface  in  the  horizontal  plane,  such  that  the  effects  of  overturning  (including  air  entrainment  and  wave¬ 
generated  turbulence)  are  at  best  very  crudely  represented.  Work  to  investigate  the  effects  of  overturning 
dynamics  on  the  statistics  is  ongoing  using  a  Smoothed  Particle  Hydrodynamics  (SPH)  model  (see  e.g. 
Dalrymple  &  Rogers,  2006). 


///i1’  ///i1’ 


Figure  2  The  non-linear  source  term 
Sni  computed  from  the  bispectrum 
obtained  from  the  observations 
(markers)  and  SWASH  model 
computations  (solid  black  line),  and 
the  linear  flux  gradient  (dashed  red 
line)  estimated  from  model  results. 

The  lower  horizontal  axis  indicates  the 
frequency,  whereas  the  upper  axis  in 
each  panel  indicates  the  normalized 
(by  the  peak  frequency)  frequency. 
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Wave  inhomogeneity  in  stochastic  models 

Operational  wave  models  are  invariably  based  on  some  form  of  the  radiative  transfer  equation  (RTE)  with 
source  terms  added  to  account  for  non-conservative  and  nonlinear  effects.  This  implies  that  the  wave  field 
is  assumed  near-homogeneous  and  near-Gaussian  so  that  individual  spectral  components  can  be  treated  as 
independent.  Although  generally  a  reasonable  assumption  for  open  ocean  applications,  over  the 
continental  shelves  and  in  coastal  waters,  inhomogeneous  and  nonlinear  effects  can  strongly  affect  wave 
evolution.  Both  nonlinearity  and  inhomogeneity  result  in  the  development  of  cross-correlations  between 
wave  components.  These  cross-correlations  contribute  to  the  wave  field  statistics,  but  can  not  be 
transported  by  the  RTE  (which  only  transports  the  energy-carrying  variance  contributions).  To  obtain  a 
more  general  model  for  the  transport  of  wave  correlators,  we  developed  a  quasi-coherent  theory  for 
stochastic  waves  (see  Smit  &  Janssen,  2013).  This  is  a  new  modeling  paradigm  to  incorporate  wave  cross¬ 
correlators  in  stochastic  models.  We  have  tested  the  new  theory  and  implemented  it  in  a  research 
stochastic  model  suitable  for  coastal  scales.  This  is  an  essential  step  toward  developing  operational 
modeling  capability  for  inhomogeneous  and  non-Gaussian  effects  near  the  coast.  Work  toward  a 
nonlinear  extension  of  the  model  is  ongoing. 

Operational  model  performance  tested  with  drifter  observations 

We  have  used  recent  observations  of  wave-current  interaction  in  the  mouth  of  the  Columbia  river  to 
validate  the  performance  of  operational  wave  models  (SWAN)  in  such  regions. 
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RESULTS 


Nonlinear  statistics  of  overturning  waves 

To  develop  a  benchmark  to  test  stochastic  closure  models,  we  have  tested  the  performance  of  a  non¬ 
hydrostatic  model  for  simulation  of  nonlinear  wave  propagation  through  a  dissipative  surfzone  (see  Smit 
et  al.  2014a).  Non-hydrostatic  models  are  very  efficient  but  do  not  capture  overturning  dynamics  in 
plunging  breakers.  To  develop  a  deterministic  modeling  platform  for  closure  validation,  which  includes 
overturning  waves,  we  have  implemented  a  Smoothed-Particle  Hydrodynamics  model  (see  Dalrymple  & 
Rogers,  2006).  The  deterministic  comparison  to  laboratory  video  observations  of  waves  breaking  over  an 
artificial  reef  is  excellent  (see  Figure  3),  suggesting  that  both  nonlinearity  and  dissipation  are  accurately 
represented  in  the  model.  Longer  simulations  and  Monte-Carlo  simulations  to  obtain  reliable  statistics  of 
laboratory  observations  (Boers,  1996)  are  ongoing. 


Observations 


SPH 


Figure  3  Smoothed-Particle 
Hydrodynamics  (SPH)  simulations 
of  waves  breaking  over  an  artificial 
reef  in  the  laboratory  (see 
Henrique2,  2004).  Top  panel 
shows  three  snapshots  of  the 
overturning  (and  breaking)  wave 
over  the  reef.  Bottom  panel  shows 
a  rendering  of  the  water  surface  as 
predicted  by  the  SPH  model  (see 
Dalrymple  &  Rogers,  2006).  The 
agreement  in  the  breaker  dynamics 
predicted  by  the  model  and  seen  in 
the  observations  is  excellent. 


Field  data  dissemination 

The  new  datasets  made  available  through  www.waveserver.org  include  data  from  the  2009-2010  High- 
Resolution  Air-Sea  interaction  DRI,  and  the  1999  ONR-funded  Shoaling  Waves  Experiment  (SHOWEX). 


Statistics  of  coastal  wave  interference 

The  quasi-coherent  (QC)  theory  developed  by  Smit  &  Janssen  (2013)  provides  a  stochastic  modeling 
platform  to  model  coherent  interferences  due  to  refraction  near  the  coast.  Using  the  coherent  length  scale 
inherent  to  the  wave  field,  we  derived  a  consistent  model  that  can  be  written  in  an  RTE-like  form  as  (see 
Smit  et  al.  2014b) 

4£  +  cx-VJT  +  ck-ViT  =  ‘V  (2) 


where  5,  denotes  partial  differentiation  with  respect  to  the  subscripted  variable,  Vx  =  ^3^,3  1  ar|d 


Vk  = 


dkt’dK 


The  cx  and  ck  represent  transport  velocities  through  physical  and  spectral  space 


respectively.  The  left  side  of  equation  (1.1)  is  an  RTE-like  transport  equation  for  the  coupled-mode 
spectrum  E  ,  and  the  right  side  of  (2)  represents  a  source  term  to  incorporate  the  development  of  wave 
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cross-correlations  (inhomogeneities)  due  to  medium  variations  (for  details  see  Smit  et  al.  2014b).  By 
recasting  the  QC  theory  from  Smit  &  Janssen  (2013)  in  the  form  of  (1),  the  similarities  (and  differences) 
with  standard  RTE-type  models  are  more  clear  and  numerical  implementation  of  QC  theory  is  simplified. 
Comparison  of  QC  and  RTE  model  simulations  to  observations  of  wave  evolution  during  NCEX,  shows 
that  in  the  presence  of  strong  refraction  and  in  regions  of  coherent  interference,  QC  theory  provides  a 
better  approximation  of  the  bulk  wave  statistics  (see  Figure  4).  In  addition  to  improved  bulk  statistics  in 
regions  where  coherent  interference  is  important,  models  based  on  QC  theory  can  provide  more  complete 
statistical  information.  For  instance,  in  regions  of  coherent  interference,  partial  standing  wave  motion 
occurs  in  the  lateral  direction  (approximately  perpendicular  to  the  mean  wave  direction  of  propagation), 
which  affects  circulation  and  transport  dynamics.  Such  information  is  not  available  in  RTE  models  (since 
wave  components  are  assumed  independent  from  the  outset).  In  contrast,  a  model  based  on  QC  theory  can 
produce  the  complete  spatial  covariance  and  thus  describe  the  complete  second-order  statistics. 


y  (km)  y  (km) 


Figure  4  Comparison  QC  and  RTE  model  simulations  to 
observations  of  wave  evolution  across  a  nearshore  canyon 
as  observed  during  the  ONR  Nearshore  Canyon 
Experiment  (NCEX),  in  2003  (see 
http:/  /  waveserver.org/ nearshore-canyon-experiment- 

ncex/).  The  comparison  is  for  a  south  swell  field  with  a 
peak  period  of  approximately  18s,  as  observed  on  16 
November  2003  (see  Smit  et  al,  2014b).  Top  panels  show 
normalized  wave  heights  for  QC  approximation  (left  panel) 
and  RTE  (right  panel).  Bottom  panels  show  normalized 
wave  height  along  15m  (left)  and  10m  (right)  depth 
contours.  In  most  regions  the  agreement  between 
observations  and  models  is  quite  good  and  comparable. 
However,  near  the  canyon  heads,  where  coherent  effects 
are  important,  the  wave  heights  predicted  by  the  QC  model 
are  in  much  better  agreement  with  the  observations  (where 
available)  than  the  RTE  model  predictions. 


An  example  of  the  spatial  covariance  function  for  the  same  wave  field  as  depicted  in  Figure  4  is  shown  in 
Figure  5,  for  two  different  center  locations  (left  and  right  panel,  respectively).  This  analysis  shows  that 
coherent  effects  are  important  near  the  Scripps  canyon  head  where  the  wave  covariance  function  develops 
nodal  lines  roughly  parallel  to  the  mean  direction  of  wave  propagation,  indicative  of  partially  standing 
wave  motion  in  the  lateral  direction  (see  also  Smit  et  al.,  2014b).  The  presence  of  standing  wave  motion 
was  confirmed  from  a  cross-correlation  analysis  of  the  available  field  observations  (not  shown  here,  see 
Smit  et  al,  2014b).  The  spatial  covariance  function  also  shows  that  just  north  and  south  of  the  Scripps 
Canyon  the  wave  field  remains  some  correlation  despite  the  complicated  refraction  and  scattering  of  wave 
energy  across  the  complex  coastal  bathymetry. 
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Figure  5  Spatial  covariance  function  of  the  same 
wave  field  as  considered  in  Figure  4.  Solid  grey  lines 
represent  wave  rays  computed  for  this  wave  field. 
Left  panel  shows  covariance  function  centered  just 
north  of  the  Scripps  canyon  head  (center  location 
indicated  by  white  cross).  Right  panel  shows 
covariance  function  centered  roughly  midway 
between  La  Jolla  and  Scripps  canyon  (white  cross). 
From  the  covariance  function  centered  north  of  the 
canyon  head  (left  panel)  we  see  that  a  nodal  pattern 
spreads  out  to  the  north  and  some  correlation  exists 
between  waves  north  and  south  of  the  canyon.  The 
covariance  function  centered  between  the  canyon 
heads  (tight  panel)  is  indicative  of  a  fairly 
homogeneous  wave  field  (for  which  the  RTE  would 
be  a  reasonable  approximation). 


The  covariance  function  centered  in  between  the  Scripps  canyon  in  the  north  and  La  Jolla  canyon  in  the 
south  is  quite  similar  to  homogeneous  theory  (not  shown)  suggesting  that  coherent  effects  are  less 
important  in  this  region.  The  enhanced  capability  of  QC  theory  to  resolve  the  covariance  function  and 
identify  standing  wave  motion  and  inhomogeneous  effects  can  help  improve  coastal  circulation  and 
sediment  transport  models,  which  rely  heavily  on  a  correct  representation  of  wave  forcing  of  the  mean 
flow  and  are  sensitive  to  spatial  variations  in  wave  height,  especially  if  the  conditions  are  close  to  the 
threshold  for  sediment  motion. 
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Impact/Implications 

The  model  improvements  developed  and  tested  in  this  study  will  contribute  to  improvements  in  modeling 
capability  of  nearshore  wave  propagation  in  research  and  operational  models.  Efficient  and  accurate 
approximations  for  four-wave  interactions  will  improve  prediction  of  spectral  shapes  in  operational 
models.  The  development  of  general  evolution  equations  for  wave  correlators  form  a  basis  for  the 
development  of  a  new  class  of  stochastic  models  that  include  inhomogeneous  and  non-Gaussian  statistics. 
In  turn,  improved  modeling  capability  of  wave  dissipation,  spectral  evolution,  and  higher-order  statistics 
such  as  skewness  and  asymmetry,  will  contribute  to  improvements  in  research  and  modeling  of  coastal 
circulation  and  transport  processes. 

Related  Projects 

The  development  of  transport  equations  for  cross-correlations  in  random  waves  also  contributes  to  the 
study  of  coastal  wave-current  interaction  as  part  of  the  ONR  Inlets  and  River  Mouths  DRI. 
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